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Abstract 

The transformation properties of strong penguin operators under the action of the flavor group 
change when they are considered as operators in (partially) quenched QCD instead of the un- 
quenched theory. As a result, additional operators and new low-energy constants appear in the 
effective theory describing non-leptonic kaon decay matrix elements in the partially quenched set- 
ting. These new low-energy constants do not have a counterpart in the unquenched theory, and 
should thus be considered as an artifact of the quenched approximation. Here we consider strong 
penguin operators consisting of products of two left-handed flavor currents, and give a complete 
one-loop analysis in the effective theory for K to vacuum and K + to n + matrix elements. We 
find that the new low-energy constants already appear in these matrix elements at leading order. 
This implies that (partially) quenched lattice computations of for instance the AI =1/2 rule are 
affected by ambiguities intrinsic to the use of the quenched approximation at leading order. The 
only exception is the partially quenched case with three light sea quarks, consistent with general 
expectations. Our results are also relevant when the charm quark is kept in the theory. 
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1 Introduction 



Recently, there has been a renewed interest in lattice computations of weak matrix elements 
relevant for the understanding of aspects of non-leptonic kaon decays, such as the AI =1/2 
rule and e'/e, which parametrizes direct CP violation in K — > ttti decays - for a recent 
review, see Ref. 0. However, while this renewed interest is due to the advent of lattice 
fermions with very good chiral symmetry, to date all such computations have been done 
in the quenched approximation, in which sea-quark effects are ignored - recent quenched 
results for both the AI = 1/2 rule and e'/e can be found in Refs. OIHj). 

In the theory where the charm quark has been integrated out penguin operators play 
an important role. In particular, referring to a commonly used basis, the left-left (LL) 
penguin operator (9 2 plays an important role in the AI = 1/2 rule, while the left-right 
(LR) operator C 6 gives a major contribution to e'/e 

When one makes the transition from unquenched QCD to partially-quenched (PQ) QCD, 
the theory is changed from the physical theory with three light quarks to a theory with 
K light valence quarks and N light sea quarks. Fully quenched QCD is the special 
case with iV = 0. This implies that the flavor symmetry group changes from the usual 
SU(3) L x SU(3) R to the graded group SU(K + N\K) L x SU(K + N\K) R 0. In general, 
this implies that the classification of weak operators with respect to the flavor symmetry 
group also changes. In particular, what happens for strong penguins is that the penguin 
operator which transformed as a component of one irreducible representation (irrep) of 
SU(3)l x SU(3) R (the octet representation) now splits into several parts, each of which 
transforms in a different representation of the PQ symmetry group. One of those is the 
"natural" generalization of the original penguin operator to the PQ theory, whereas the 
other transforms in a more complicated way under SU(K + N\K)i x SU (K + N\K) R . We 
will refer to these two parts as the "singlet" and "adjoint" operators, respectively - for 
reasons that will become clear in Sec. 2. 

At the level of the effective theory, this means that new low-energy constants (LECs) 
occur for the adjoint operators, with no counterpart in the unquenched theory, and which 
thus should be considered an artifact of the use of an approximation where the number 
of light sea quarks is not equal to three. 1 These new LECs appear in physical matrix 
elements, unless one decides to drop the corresponding adjoint operators from the PQ 
theory - meaning that the definition of the penguin operators themselves is changed in the 
transition to the PQ theory. We observe that also the LECs for the singlet operators do 
not have to be equal to those of the physical three-flavor theory if iV ^ 3; even their scale 
dependence will in general be different. 

This problem was considered in a previous paper for the strong LR operators 0$$ jB] . Here 
we consider the same problem for the LL operators 0\^- While the group theory involved in 
the LL case is a little more complicated than for the LR case, our results are rather similar. 

■""We will see that for N = 3, and when choosing valence- and sea-quark masses equal, physical matrix 
elements of all new operators vanish. 
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We find that also in the LL case the new adjoint operators do contribute to physical matrix 
elements already at the leading order (at tree level) in chiral perturbation theory (ChPT). 2 
Here we demonstrate this explicitly with the examples of the K° to vacuum and K + to 
7r + matrix elements. We also find that the one-loop corrections for the adjoint operators 
differ from those of the singlet operators, calculated in Ref. J2]. In other words, the singlet 
and adjoint LECs do not occur in some fixed, given linear combinations in physical matrix 
elements beyond tree level. 

The outline of this paper is as follows. In Sec. 2 we show how the standard LL strong 
penguins break up into terms that transform differently under the enlarged symmetry 
group of PQ QCD with K valence quarks and N sea quarks. In Sec. 3 we construct the 
corresponding weak operators in ChPT, while Sec. 4 contains a discussion of some relevant 
peculiarities that arise for representations of the graded symmetry group SU(K + N\K) 
in the PQ theory. Section 5 contains explicit expressions for the K° to vacuum and K + 
to 7r + matrix elements to one loop in ChPT. Section 6 summarizes our results for the 
fully quenched case, and the final section contains our conclusions. Some group-theoretical 
details and useful relations amongst weak operators are given in two appendices. 



2 Left-Left penguins in partially quenched QCD 

We consider the LL penguin operators 

01 = (sd) L (uu) L - (su) L (ud) L (2.1) 
= (s a d a ) L (upUp + dpdp + spSp) L - {s a dp) L (upu a + d a d p + s p s a ) L , 

2 = (sd) L (uu) L + (su) L (ud) L + 2(sd) L (dd + ss) L 

= (s a d a )L(Uf3U{3 + dpdp + SpSp) L + (Sadp) L (ilf3U a + d a d(3 + ~SpS a ) L , 

where 

(QiQ2)l = lilnPLq2 , (2.2) 

with the projection operator on left-handed chirality Pl = (1 — 7s)/2. In the second expres- 
sion for each operator, we have made the color indices a,/3 explicit. Both operators 0^2 
are penguin operators, and each is a linear combination of color un-mixed and color mixed 
terms. Both operators transform in the octet representation of SU(3)l, and, trivially, in 
the singlet representation of SU(3)r. They are part of a basis of irreducible representa- 
tions of the chiral group that are CPS invariant and with definite isospin 7=1/2 and 
1 = 3/2 jlj. In Appendix A we clarify how this basis is related to a set of weak operators 
more frequently used in phenomenological analyses. Our basis is especially convenient for 
working out group theoretical properties. 

As already mentioned in Sec. 1, when we consider the LL penguin operators 0^2 of Eq. (J2.1)) 
in the partially quenched theory, the representation content of these operators changes. A 
general realization of PQ QCD contains K valence quarks, each accompanied by one of 

2 A speculation in Ref. 6 that they would only start contributing at order p 4 turns out to be incorrect. 
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K ghost quarks with the same mass in order to suppress the valence- fermion determinant, 
and TV sea quarks - the dynamical quarks - which can all have masses different from those 
of the valence quarks. The relevant flavor symmetry group enlarges from the physical 
SU(3) L x SU(3) R to the graded group SU(K + N\K) L x SU(K + N\K) R 0. 3 

It is clear that the (sg?)l factors of both operators in Eq. (|2.1j) are still a component of 
the adjoint representation of SU(K + N\K), 4 while the factors (uu + dd + ss)l no longer 
transform as singlets. Instead, the operators can now be written as 

0l = K PQS + PQA^ (2-3) 



N 

Q 2 = K P QS + PQA 



(D± QS = (q a Aq a ) L {q^)L ± {q a ^qp)L{q p q a )L , 

0± QA = (q a kq a )L{qpMp)L ± (q a kqp)L{qpM«)L , 

where we introduced the spurion fields A and A with values 

V = <W* i2 > (2.4) 
( K K 

= diag { 1 -N---~N- 

The first K diagonal elements of A are equal to (1—K/N) - corresponding to the K valence 
quarks - and the last N + K diagonal elements are equal to (-K/N) - corresponding to 
the iV sea quarks and the K ghost quarks, both of which do not occur in Oi >2 - The quark 
fields are graded vectors in flavor space, with fermionic components given by the valence 
and sea quarks, and bosonic components by ghost quarks. The indices i and j are graded 
flavor indices, and run over valence, sea and ghost flavors. 5 For the down (strange) quark 
we have i — 2 (i — 3). The completely quenched theory, i.e. the theory with N = 0, for 
which the split of Eqs. (|2.3U2.4j) is singular, will be dealt with in Sec. 6. 

We note at this point that for iV = K = 3 we regain the physical three-flavor theory. 
The adjoint operators Q^ A now contain only terms involving either sea or ghost quarks, 
and it is rather straightforward to see that their contributions to physical matrix elements 
{i.e. those with only valence quarks on the external lines) vanishes because of cancellation 
between sea-quark and ghost-quark loops. For this cancellation to happen, valence masses 
and sea masses should be chosen equal. 

The spurions A and A both transform in the adjoint representation of SU(K + N\K)l, as 
can be seen from the fact that both have a vanishing supertrace (str) JU]- The operators 

3 For a detailed analysis of the actual symmetry group in the euclidean lattice theory, we refer to 
Ref. The upshot is that for our purposes, it is appropriate to consider the PQ symmetry group to be 
SU(K + N\K) l x SU(K + N\K) R 0. 

4 We will often drop the subscript L on the group from now on, since all operators considered in this 
paper are trivial with respect to the right-handed group. 

5 Wherever explicitly written, we denote flavor indices with color indices with a, (3, . . . and Dirac 

indices with a,b, . . .. 
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thus transform in the adjoint representation, while the operators Q±® A transform 
as the product representation of two adjoint irreps, and they are thus reducible. The cor- 
responding decomposition of Q^ A is accomplished by (anti-)symmetrization in covariant 
and contravariant indices, and by "removing" supertraces on pairs of covariant and con- 
travariant indices, much as is done in the case of SU(N) ^Oj- Here we take the quark 
fields qi as covariant, and the anti-quark fields <f as contravariant. It turns out that the 
operators Q^ A and 0^® s (0+® A and 0+® s ) are already symmetric (anti-symmetric) in 
both their two covariant and their two contravariant flavor indices - see Appendix B for 
details. 

In the next section, we are going to construct the low-energy bosonized effective lagrangian 
for the operators in the PQ theory. For this we do not need the details of the decom- 
position into irreps, and we postpone further discussion of this decomposition to Sec. 4. 



3 The effective lagrangian 

The bosonized low-energy effective lagrangian is constructed in terms of the non-linear 
field 

£ = exp(2i<V/) , (3.1) 

where $ is the (2K + N) x (2K + N) hermitian field describing the Goldstone mesons, and 
/ is the pion-decay constant in the chiral limit (normalized such that f n = 132 MeV). Out 
of this field, we construct the necessary building blocks 

Lfj, = ffifyEt , (3.2) 
X+ = 2fi (SM t + MS f ) , 

where M is the (2K + N) x (2K + N) quark-mass matrix, and 2B = —(uu)/f 2 in the 
chiral limit ^Tj. These building blocks as well as the spurions A and A all transform 
in the same way under the left-handed group SU(K + N\K), and to lowest order in the 
chiral expansion, we can construct the effective operators (note that [A, A] = and that 6 
str (L M ) = 0) 

£ A = str (AL M ) str (AL M ) , (3.3) 
C A = str (AL^AL^) , 
C A = str (AAL^) , 
C A = str (AAX + ) . 

In deriving this list, we have also used CPS symmetry [12J. For the singlet operators 0±® s 
one replaces A—*l, and the above operators reduce to 

Cf = str (AL M L M ) , (3.4) 
£$ = str (AX+) . 



3 If we consider the PQ theory with the 77' integrated out. 
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The singlet operators Q±® s are represented in chiral perturbation theory by the leading 
order lagrangian O E] 

C s = -aW £f + 4 8 ' 1} C* , (3.5) 

where are two weak LECs. 7 Note that the LECs for the two different operators 0±® s 
are independent of each other, even though we use the same symbol for both of them. 

The operators Q^ A correspond to different representations of SU(K+N\K), and are thus 
represented by different linear combinations of £1234- A mostly straightforward analysis 
leads to the bosonization rules 8 

O p ± QA ^C A = a A ±(£ A ± C A ) + a A ±£ A + a A± C A . (3.6) 

Here we have explicitly indicated the dependence of the LECs on the operator through the 
superscripts ±, because they refer to different representations of the PQ flavor group. We 
conclude that the transition from the unquenched theory to the PQ theory leads to the 
introduction of three new LECs for each of the two operators 0\ and Oi- 



4 Mixing of four- quark operators 

As stated already in the previous section, the two operators Q±® A correspond to two 
different representations of SU(K + N\K). Both four-quark operators may be written in 
the form 'q^q'T^'q^qi, with T symmetric (or anti-symmetric) in both i j and k <-> /. 9 
Further decomposition of these representations is generally possible by splitting the tensor 
T into a part which is supertraceless on any pair of covariant and contravariant indices, and 
a supertrace part. Note that, because of the value of the spurion A, the double supertrace 
of T on both pairs of indices vanishes in our case. 

Taking the supertrace of T on one pair, one obtains a tensor , and one may thus construct 
new four-fermion operators by replacing the tensor T by the new tensor S^Sj, of course 
after (anti-) symmetrizing this new tensor in correspondence with the symmetry properties 
of the tensor T. This results in the two new four-quark operators, 

Ol QT = (q a AAq a ) L (qpqp) L ±(q a AAq (3 ) L (q (3 q a ) L (4.1) 

Our task is now to see whether these new operators are contained in the original four-quark 
operators (D±® , i.e. whether the operators 0±^ A and 0^ T mix. We observe that these 

7 The minus sign is there to make the definition of ajj 8 ' 1 ^ conform to that of Ref. 12 , in which the 
effective lagrangian was defined in Minkowski space. We work in euclidean space. 

8 Some of the less straightforward aspects will be discussed in the next section. See also Appendix B. 

9 (Anti-)symmetrization here is understood as appropriate for representations of graded groups. "Sym- 
metrization" implies symmetrization in bosonic indices and anti-symmetrization in fermionic indices, and 
vice versa (see Appendix B and Ref. [TU|). 
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operators look just like the singlet operators 0± QS , but for the factor 1 — K/N, which is 
the only remnant of the spurion A. 

There are three types of diagrams when one considers QCD corrections to all these op- 
erators: self-energy corrections on the external lines, vertex corrections, and penguin-like 
diagrams. The mixing is only possible through the penguin-like diagrams, because the first 
two types of diagrams do not change the flavor structure of the operator. This means that 
we have to consider what happens when we contract one quark field with one anti-quark 
field in each of these operators. Such a contraction, in group-theoretical language, corre- 
sponds to taking a supertrace on a pair of quark and anti-quark indices of the tensor T. 
Performing the contraction inside the operators 0±® A , we obtain 10 




where S^g is the quark propagator for flavor q in an arbitrary gluon background - i.e. the 
quark propagator with an arbitrary number of gluons attached - and we used that 

ltl P L Sl pl ,P L = -^P L tr (^PlS^) . (4.3) 

The latter relation is easily proved by expanding the left-hand side on a basis of 2 4 euclidean 
hermitian Dirac matrices for fixed color indices a and (3. In these equations, tr stands for 
a trace over Dirac indices only. 

Quark masses in our theory can be thought of as insertions leading to higher- dimension 
weak operators, and we may thus consider the same result ()4.2|) in the chiral limit, in 
which case the fermion propagator S q — > S becomes flavor independent. Equation ()4.2|) 
thus simplifies to 

PQA _ _ 2 A ^ ((s a ^P L d ) tr (^P L S a p) ± (Sa^P L d a ) tl (^P L S^)) . (4.4) 

The above expression vanishes if we take K = N, consistent with the fact that the operators 
qPQT vanigh j n this case (c/. Eq. (|4.1j) ). When K ^ N, this result shows that indeed the 
operators 0± mix with O p ± qt , because the gluons attached to the fermion propagator 
couple to the flavor-singlet bilinear g7 M P L g. 

10 Wick contraction of graded quark bilinears gives q^qia = (—l) 9 ^ S^Sj, with index g(j) denned in 
Appendix B. 
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Performing the same contraction on the operators 0± , we find 

O pqt ^ (l " f ) (TN - 2) {{s alil P L dp) tr { ltl P L S aP ) ± (s al ,P L d a ) tr (7^^)) . 

(4.5) 

Comparison of Eqs. (|4.5jl and (|4.4jl implies that there exist the supertraceless linear com- 
binations 

that do not mix with the operators (D± . In group-theoretical language, these linear 
combinations are irreducible under SU(K + N\K). We thus find that in general, except 
for K = N, our four-quark operators do "contain" the operators V± . In the effective 
theory, the latter correspond to the operator Cf in Eq. ()3.3j) . justifying the bosonization 

A closer look at Eq. (|4.fi|) reveals that when the number of sea quarks, N, is equal to two, 
a singularity arises when one tries to decompose O p ® a into irreps. Clearly, for N = 2 
O pqa does mix with O p ® t , but it is not possible to define an operator which does not 
mix with P _ QT . In group theoretical language, the corresponding statement is that the 
representation in which (J P ® A transforms is not fully reducible. This is indeed what one 
finds if one analyzes the tensor Tjj for this operator, and in accordance with the fact that 
such reducible but not decomposable reprentations do exist for graded groups [EH- 11 In any 
case, we still have to include the operator C A in the effective theory, and the bosonization 
rule ()3.6|) thus stays the same. 



5 K° —> vacuum and K + — > ir + matrix elements 



In this section we give results for the simplest kaon matrix elements of the new weak 
effective operators, £± (cf. Eq. ()3.fij) ). We first give the tree- level results, and then include 
also the chiral logarithms which occur at 0(p 4 ). We will not present a detailed analysis 
of 0(p 4 ) contact terms, because it is unlikely that the matrix elements of O p ® a will be 

(8 1) 

numerically computed in the future. The reason is that only the singlet LECs a\ 2 are 
the interesting ones for physical predictions, and in a PQ setting they can be obtained 
from the operators 0±® s , as long as the number of light sea quarks is physical, i.e. N = 3 
P^j . For a complete 0{p 4 ) analysis of the singlet operators represented at lowest 

order by C s in Eq. (|3.5p . we refer to Ref. [7] (see also Ref. [13]). On the other hand, it 
is relevant to verify whether or not the tree level and one-loop contributions to physical 
matrix elements of C A have the same form as those arising from the singlet operator C s . 

At tree level, we find that 

(0|£^> = W&ZM U _ ^ at , (5.1) 



/ V AT 



11 The operator Q P ^ A corresponds to a tensor which is symmetric in both % < — ► j and k «-> I, with the 
appropriate grading (see Appendix B). 
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where in the case of K + — > n + we took the kaon and pion masses to be equal, Mk = 
M w = M. It is clear that the new LECs oc^\ b 2 already appear at leading order in ChPT, 

competing with the leading-order contributions coming from the singlet LECs a\ ' 2 . 

Next, we give the non-analytic terms arising at 0(p 4 ). We express our results in terms of 
bare meson masses, where we take all sea quarks degenerate for simplicity, and also work 
in the isospin limit, setting m u 

(5.2) 



= m d 


= m, i.e. 


M K 


= B (m + m s ) , 


Ml 


= 2B m , 


Mis 


= B (mi + m sea 


Mis 


= 2B m seSL , 


Ml 


= B (mi + mj) 



where m ; is the mass of the z-th valence quark - we use the labels i = u,d,s respectively 
i = 1,2,3 interchangeably to label valence quark masses. For simplicity, we set Mk = 
M n = M in the K + — ► n + matrix elements, thus working in the degenerate limit. It is 
further assumed that the rf of the PQ theory is heavy, and therefore has been integrated 
out [ini HZ1 El E| • We give MS expressions for all one-loop results. 

For the non-analytic terms in (0\O±® A \K°) at one loop we find 
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(5.3) 



where m 3 «-> m 2 stands for the exchange of strange and down valence-quark masses, and 
we took the exact isospin limit. 
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For the non-analytic terms in (ir + \0± at one loop, with Mk 

M$ s = (M 2 + M 2 SS ) /2, we find 
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(5.4) 



Eqs. ()5.3|) and (|5.4|) represent the partially quenched one-loop contamination which affects 
K to vacuum and Af + to tt + amplitudes through the operators 0±® A . By construction, 
these contributions vanish in the full QCD limit, when N = K and the valence- and sea- 
quark masses are chosen to be pairwise equal. We note that in Eq. (|5.4|) M denotes the 
mass of a meson made out of two valence quarks, i.e. M = for all i, j valence. The 
contribution of sea-quarks can be traced through the index S on M vs , which is the mass 
of a meson made out of one valence and one sea quark, and M$s, which is the mass of a 
meson made out of two sea quarks. 



6 The quenched case 

In this section we summarize the results for the quenched case, i.e. the case without any sea 
quarks, for which the symmetry group is SU(K\K). The quenched case needs to be treated 
separately, because the rf cannot be integrated out [TBI EH] • Moreover, the phenomenon 
that not all reducible representations of graded groups are fully reducible already shows 
up when one tries to construct the adjoint representation: the product of the fundamental 
irrep and its complex conjugate is not fully reducible The decomposition (J2.3j) now 
takes the form 

£>i = — *— (-y0^ s + 0^ NS ) , (6.1) 
1 + 7 v ' 

o 2 = -J- (jof + o2 NS ) , 
1 + 7 v 7 

0± S = (? a Ag Q ) L (g /3 g /3 ) L ± {q a Aqp) L (q p q a ) L , 
02 NS = {q a Aq a ) L {q p Nq p ) L ±{q a Aq p ) L {q Nq a ) L , 
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where the spurion field N is given by 

N = diag (1, . . . , -7, . . .) ■ (6.2) 

The first K diagonal elements of N are equal to 1, and the last K diagonal elements 
are equal to —7, with 7 arbitrary except 7 7^ — l. 12 Note that N has a non-vanishing 
supertrace, unlike A in the PQ case, while the supertrace of the unit matrix diag(l, . . . , 1) 
does vanish, again unlike the PQ case. The operators 0± s and 0± NS each are represented 
by their own LECs in the effective theory, just as in the case of LR penguins 6j. Naively, 
it looks like one could avoid those LECs which correspond to 0± by choosing 7 = 0. 
However, the operators 0± mix with 02 S for any choice of 7, as can be shown in the 
quark language using arguments similar to those of Sec. 4. 

Once again, in order to analyze results from quenched QCD, one has to develop the ef- 
fective theory for the operators 0± NS . This is straightforward: one exactly follows the 
construction of Sec. 3 for the PQ case, by replacing the PQ spurion A everywhere with N. 
However, there is one extra operator, because in the quenched case str 7^ - due to 
the rj', which cannot be integrated out in the quenched case - and the quenched non-singlet 
effective lagrangian replacing Eq. ()3.6|) reads now (see also Appendix B) 

Cf s = a?* (str (AL M ) str (NLJ ± str (AL M A>L M )) (6.3) 
+a$ k (str (KNL^L^) ± str (AA>L M ) str (L M )) + a^str (AA>X + ) . 



The tree-level K to vacuum and K + — > 7r + matrix elements are given by 

2 

<__ 

f 



C» ± \K°) = il{M « ~ M -> , (6.4) 



AM 2 

[it \C±\K } = (Ta la -a lb - a 2 ), 

where we labeled the non-singlet effective operators and LECs with a superscript N instead 
of a superscript A for the quenched case. 

Before presenting quenched one-loop results, we remind the reader that the tree-level prop- 
agator of a neutral valence meson made out of quark and anti-quark flavor % is given by 

£>..M = — ^ - m ° + ap2 (6 5) 

where m\ is the "double-hairpin" vertex at zero momentum, and a parametrizes its mo- 
mentum dependence |18j . For the non-analytic terms at one loop, we find 

(OlO^i^e-Zoop = ^f,K 4 3 fl-l0g#l-J(l + 7) £ 



2 In Ref. |S] we chose 7 = 1. 
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M* (ml - «M 2 ) (l - log fj - M 3 4 3 (m 2 - «M 3 2 3 ) ^1 - log 

1 /2 \ M 2 

±-M 2 3 (m 2 - 2aM 2 3 ) =F M| 3 (^-m 2 - aM 2 3 J log -f 



±(l+ 7 ) £ M 3i l-lo{ 



i valence 



- (m 3 <-> m 2 ) 



4ia 



iV± 
16 



1 

— m; 



o 



M 3 2 3 - (|m* 



±M 3 4 3 (l - log^ - (m 3 <- m 2 )| 
2/n " (M|-M 2 )jm 2 -3aM 2 3 + 



aM 2 3 M 2 3 log 



Mf 
/i 2 



3/(4tt/) 2 
+ (m 3 <-> m 2 ) ^ , 



-m 2 + 2aM 2 3 -2M 3 2 3 -£ 



M| 3 - M 9 2 < 



22 



log 



Mis 

(6.6) 



for the (OIO^^I^ ) matrix element, and 

4M 2 < ± 



(71 

—loop 



f2( j ;t . )2 ,2(6M 2 + (1 + 7 )(m 2 - 3aM 2 )) log 



M 2 

7 2 



-6M 2 + K(l + 7)(m 2 + aM 2 ) ± (2m 2 + (4a + 3X(1 + 7))M 2 ) 
± (6m 2 - (16a + 3K(1 + i))M 2 ) log ^ j 

SM 2 a^ \ , 2 2N1 M 2 2 2 2 / , M 2N 

16 J -(3m 2 - 8aM 2 ) log — - m 2 - 2aM 2 ± 3M 2 1 - 2 log — 



3/ 2 (4tt/)2 
8M 2 af 



M 2 

m 2 + (m 2 - 4aM 2 ) log — 



3/ 2 (4tt/)2 
for the (7r + |(92 Ar,S '|-ft' + ) matrix element. 



(6.7) 



7 Conclusions 

We add a few comments to the main conclusions already stated in the Introduction. 

First, a consistency check on all our explicit tree-level and one-loop results is that they 
should vanish when the number of valence quarks K is chosen equal to the number of 
sea quarks N and when the valence- and sea-quark masses are also chosen to be pairwise 
equal. We noted this already in Sec. 2, and it is straightforward to see that indeed all 
results contained in Sec. 5 do satisfy this requirement. 

We emphasize that the new adjoint operators 0± occurring in the PQ theory are gen- 
uinely new operators, and one thus expects that one-loop corrections in ChPT for matrix 
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elements of these operators differ from those of the singlet operators 0±® s . We find that 
this is indeed the case. If one however only considers the tree-level results, only one particu- 
lar linear combination of singlet and adjoint LECs appears in all matrix elements (including 
K —>■ 7T7T matrix elements). This is unlike the case of LR penguins, where enhancement of 
the adjoint operators leads to the appearance of chiral logarithms already at leading order 
in ChPT [Bj. Thus, Refs. 0111], while aiming for the leading-order LECs for 0\ and O2 
in the quenched approximation, have actually computed the linear combinations of LECs 
that appear at tree level, namely 

-L- M 8 - 1 ) + <±) , (7.1) 
1 + 7 v ' 

1 / (8,1) _,_ jV± JV± / (8,1) . N±\\ 

1 + 7 v * 

At first, the appearance of the arbitrary parameter 7 looks a bit strange. However, this 
just reflects the fact that the non-singlet operators Q% NS are only defined modulo mixing 
with the singlet operators 0± , since they do not by themselves constitute an irrep, as 
discussed in Sec. 6. In the PQ case no such arbitrariness arises. 

Our results show that quenching artifacts do modify, already at the leading chiral order, 
those AS* = 1 weak matrix elements that receive contributions from LL penguin opera- 
tors. 13 This is especially the case for the AI = 1/2 rule where the dominant contributions 
come from the current-current operators Q± and Q2 of Eq. (jA.2|) . In addition, quenching 
contaminations to LL penguin operators can in principle affect a lattice determination of 
e'/e. 14 While quenching artifacts analyzed in Ref. jE] affect a dominant contribution to 
e'/e coming from the LR penguin operator Q s , the quenching ambiguity affecting the op- 
erator Q4 through 0\ and O2 can be relevant in the presence of a large cancellation of the 
dominant contributions from Q 6 and the electroweak penguin operator Q 8 [2*T] . 

Finally, we note that our results are also relevant for the theory in which the charm quark 
is kept. In that case, the relevant weak operators are 0±, which can be written as 

0_ = {sd) L {uu) L -{su) L {ud) L -{u^c) (7.2) 

= O x - (u -> c) , 

+ = (sd) L (uu) L + (su) L (ud) L - (u -> c) 
1 2 , _ 2 , « 

= T®2 + J^0 3 + -04 - (U - C , 

5 15 3 

where 0^ ; 4 transform in the 27-dimensional irrep of SU(3)l - see also Appendix A. The 
operators appear, and the discussion of this paper thus also applies to the operators 
0±. The charm quark transforms as a singlet under both SU(3) and SU(K + N\K), and 

13 We disagree with the conclusion of Ref. that the AI — 1/2 rule is not affected by (partial) quench- 
ing artifacts at tree level. In particular, we point out that it is not possible to decide unambiguously within 
the quenched approximation what is the "best" linear combination of singlet and non-singlet operators to 
choose. Only an unquenched computation can, in hindsight, decide this issue. 

14 We thank the referee for pointing this out. 
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the "(w — > c)" terms in Eq. ()7.2j) thus transform in the adjoint representations of both 
SU(3) and SU(K + N\K). Our observations here do not apply to the SU(4) case of an 
unphysically light charm quark, because the SU(4) transformation properties of 0± are 
different. 
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Appendix A: AS = 1 weak operators 

In this appendix, we clarify the relation between the operator basis we used and a set of 
weak operators often used in phenomenological analyses of kaon decays - a comprehensive 
review can be found in Ref. [22J. Our basis of AS = 1 four-quark effective weak operators 
corresponds to irreducible representations of the chiral group SU(3)l x SU(3)r, that are 
CPS invariant J2| invariant under the product of CP and the exchange of strange 

and d quark), and with definite isospin, 1 = 1/2 and / = 3/2: 



(8l, Ir) 


I 


= 1/2 


o x = 


(sd) L (uu) L - (su) L (ud) L , 




I 


= 1/2 


o 2 = 


(sd) L {uu) L + isu) L {ud) L + 2(sd) L (dd + ss) Ll 


(27 L , Ir) 


I 


= 1/2 


Os = 


{sd) L [{uu) L + 2(dd) L - 3(ss) L ] + (su) L (ud) L , 




I 


= 3/2 


o A = 


(sd) L [(uu) L - (dd) L ) + (su) L (ud) L , 


(8l,1r) 


I 


= 1/2 


o 5 = 


(sd) L [(uu) R + (dd) R + [ss) R ] , 




I 


= 1/2 


o G = 


(~s a dp) L [(upu a ) R + {dpd a ) R + (s/3S a ) R } , 


(8l,8 r ) 


I 


= 1/2 


o 7 = 


{sd) L [(uu) R - {ss) R ) - (su) L (ud) R , 




I 


= 1/2 


Os = 


(Sadf3)L[(Uf3U a ) R - (Sf3S a ) R ] - (s a U{3) L (ilpd a ) R 


(8l,8r) 


I 


= 3/2 


o 9 = 


{sd) L [(uu) R - (dd) R ] + (su) L (ud) R , 




I -- 


= 3/2 : 


010 = 


(s a d/3) L [(upu a )R - {dpd a ) R ] + {s a up) L {upd a ) R 



(A.l) 



The operators C 5 and 6 are the LR penguin operators considered in Ref. 0. A set of 
AS = 1 four-quark effective operators frequently used in phenomenological analyses of 
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kaon decays are the Qi, i 



1, 10. They are related to our basis as follows: 



g 2 



Qi 




Qa 

g 5 



Q 3 



±(o 7 + o 9 ), g 8 = I(o 8 + o 10 ) 



Q9 




Q10 



— X - —0 2 + —0 3 + -04 , 

2 10 10 2 ' 



(A.2) 



where we have used Q 4 = + Qi> Q9 — 3/2Qi — 1/2Q 3 and Q10 = 3/2Q 2 — I/2Q4 = 
Q2 — 1/2Q 3 + 1/2 Qi- It appears that any modification induced by partial quenching of 
the LL penguin operators 01,2 considered in this paper will affect the operators Qi,2, £?3,4 
and Qg,io according to the decomposition in Eq. (jA.2|) . The change of the current-current 
operators Qi^ is relevant to any (partially) quenched lattice calculation of the AI = 1/2 
rule. A change of the LR penguin operator Q 6 can drastically effect a (partially) quenched 
lattice calculation of e'/e as has been analyzed in Ref. j^j. The fact that Q4 is affected as 
well may also be relevant for the determination of e'/e. 

Appendix B: Graded groups 

In this appendix, we give a few details of graded-group technology as applied to the PQ 
operators considered in this paper. 

First, our definition of the supertrace of a matrix is 



where g(i) = 1 if the index i is fermionic (corresponding to a valence or sea quark), and 
g(i) = if the index i is bosonic (corresponding to a ghost quark). The extra overall minus 
sign is there to make the supertrace reduce to the normal trace in flavor space in the case 
of unquenched QCD. 

Next, we derive in some detail our claim that the operators 0±® s and are already 

symmetric or antisymmetric in the pairs of flavor indices of the quarks (which we refer to 
as covariant, following Ref. JU]) and anti-quarks (which we refer to as contravariant). 



str(M) = -^(-l)^M/, 



(B.l) 
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Begin with considering the operator 

{laaM'aaWpbM'pb) = (s a ad' aa )(q pb Aq' pb ) (B.2) 

in which a, (3 are color indices and a, b are Dirac indices with summation convention for 
all explicit indices, and we abbreviate 

q' = l„P L q (B.3) 
for fixed fi. Symmetrizing this in both covariant and contravariant flavor indices leads to 

{qaa^fl'aaWpbM'pb) + (9aa A ^J&) fafilMaa) + fajllMaa) {laa^b) + pb A q'p b ) (q aa Aq' aa ) 

= Hq aa ^<laa)(qpbM'f3b) + ^aM'pbWpbM'aa) ■ (B.4) 

We used that the symmetric product of two quark fields is given by ^U] 

quq2j + q2iqij = quq-ij + (-l) g{i)gij) qijqii , (B.5) 

and likewise for anti-quark fields and anti-symmetric products. We now wish to write this 
in a form in which we can use the shorthand (|2.2|) . We rewrite the second term on the 
right-hand side of Eq. (jB.4|) as 

{q aa hq'p h ){q pb Aq' aa ) = ^(-l)^^^)^/ (B.6) 

i 

= -(q a Aqp) L (qpAq a ) L , 

where in the last step we fierzed the operator, taking into account that this involves an 
extra minus sign if q ai and are bosonic {i.e. ghosts). We conclude that symmetrizing 
the operator of Eq. ()B.2j) in both quark and anti-quark flavor indices yields the operator 

qPQA 

A similar argument shows that anti-symmetrizing in both quark and anti-quark flavor 
indices leads to the operator 0+® A , and that operators with mixed symmetry (symmetric 
in quark flavor indices and anti-symmetric in anti-quark flavor indices, or vice versa) vanish. 
The operators 0±^ s and 0±® T are obtained in the same way. Note that these last two 
pairs of operators are in principle different, even though for the particular value of the 
spurion field A relevant for this paper (cf. Eq. O) we have that 0± QT = (l - f ) <D± QS . 

We may now bosonize these four-quark operators following standard techniques, see for 
example Ref. j3|. The only difference with respect to the usual, non-graded case is that 
care has to be taken with extra signs due to the grading of our symmetry group, as discussed 
briefly above, and in much more detail in Ref. pUj- This leads to the weak lagrangians 
given in Eqs. ()3.5j) and (|3.6jl . In particular, bosonization of 0±® T leads to In fact, the 
effective operators corresponding to C£ QT are 

str (AAL^L^) ± str (AAL^) str (L M ) , (B.7) 

but we have that str (L M ) = in the PQ theory with the rj' integrated out, leading to Cf 
in Eq. (J3.6|) . In the quenched theory the spurion A is replaced by N, and the 7/ cannot be 
integrated out, whence £3 of Eq. (|3.6|) is replaced by the quenched version 

str (ANL^) ± str (ANL^) str . (B.8) 
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